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An overview of the hydrodynamic profiles derived by kinetic theory tools 
(Boltzmann equation and kinetic models) for the gravity-driven Poiseuille 
flow with particles colliding either elastically (planar and cylindrical geome- 
tries) or inelastically (planar geometry) is presented. The results obtained 
through second order in the gravity acceleration show that the Navier-Stokes 
predictions are qualitatively incorrect over distances of the order of the mean 
free path, especially in the case of the temperature profile. 



1 Introduction 

In the last decade, the Poiseuille flow generated by the action of a uniform 
longitudinal body force mg (e.g., gravity), rather than by a longitudinal pres- 
sure gradient, has received much attention from computational [1, 2, 3] and 
theoretical [3, 4, 5, 6, 7, 8, 9, 10, 11, 12] points of view. This interest has 
been fueled by the fact that the gravity-driven Poiseuille flow provides a nice 
example illustrating the limitations of the classical Navier-Stokes (NS) de- 
scription, even in the bulk domain (i.e., far away from the boundary layers), 
over distances of the order of the mean free path. In contrast to what hap- 
pens with the structure of a plane shock wave, where the NS discrepancies 
are essentially quantitative and are widely accounted for by the Burnett de- 
scription [13], here the main failures of the NS predictions are qualitative and 
remain in the Burnett theory. 

Kinetic theory analyses of the gravity-driven plane Poiseuille flow based on 
an expansion in powers of the gravity acceleration g [5, 6], on Grad's moment 
method [3, 8], or on an expansion in powers of the Knudsen number [10], 
show that to second order in g the temperature proflle includes a positive 
quadratic term in addition to the negative quartic term predicted by the NS 
(and Burnett) hydrodynamic equations. As a consequence of this extra term, 
the temperature does not present a flat maximum at the middle of the channel 
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but instead exhibits a bimodal shape with a local minimum surrounded by two 
symmetric maxima at a distance of a few mean free paths. The Fourier law 
is dramatically violated since in the slab enclosed by the two maxima the 
transverse component of the heat flux is parallel (rather than anti-parallel) to 
the thermal gradient. The kinetic theory prediction of a bimodal temperature 
profile in a slab has been confirmed by computer simulations [2, 3, 7]. A similar 
behavior occurs in the case of the Poiseuille flow in a pipe [9, 11]. 

For conventional gases under terrestrial conditions the interest of the 
Poiseuille flow induced by gravity is rather academic. However, this is not 
necessarily so when dealing with a "granular" gas [14], i.e., a collection of a 
large number of discrete solid particles (or grains) in a fluidized state such 
that each particle moves freely and independently of the rest, except for the 
occurrence of inelastic binary collisions. A kinetic theory study of a gas of 
inelastic hard spheres in a slab under the action of a longitudinal force and 
excited by a white noise "heating" has recently been undertaken [12]. 

The aim of this paper is to provide an overview of the hydrodynamic pro- 
files derived by kinetic theory tools (Boltzmann equation and kinetic models) 
for the gravity-driven Poiseuille fiow in the cases of particles colliding either 
elastically (planar and cylindrical geometries) or inelastically (planar geome- 
try). The problem is presented in Sect. 2 and the NS solution to order is 
given in Sect. 3. The kinetic theory description for elastic and inelastic colli- 
sions is provided in Sects. 4 and 5, respectively. Finally, the main conclusions 
are briefly presented in Sect. 6. 

2 Statement of the Problem 

Let us consider a monodisperse dilute gas of identical spherical particles of 
mass m. The interaction among particles is assumed to reduce to binary col- 
lisions which can be either elastic or inelastic. In the first case, the (local) 
kinetic energy per unit volume is preserved by collisions. On the other hand, 
if the collisions are inelastic, energy is dissipated into the internal degrees of 
freedom; we will assume that this internal energy sink is counterbalanced by 
an external energy source produced by some sort of driving (e.g., boundary 
vibrations, external thermostats, ...). In addition, the system is under the 
influence of a constant gravitational fleld characterized by the acceleration 
vector g. Under the above conditions, the macroscopic balance equations for 
the local densities of mass, momentum, and energy read 



Dtn + nV • u 



0, Au + 




(1) 



AT + |- (V • q + P : Vu) = -(C - 7)r . 



(2) 
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Fig. 1. Sketch of the longitudinal section of the Poiseuille flow induced by a gravi- 
tational force. In the planar geometry, the boundaries are infinite parallel plates and 
the variable s represents the Cartesian coordinate y. In the cylindrical geometry, the 
boundary is an infinite cylindrical surface and the variable s represents the radial 
coordinate r 



In these equations, Di = i9t + u • V is the material time derivative, n is 
the number density, u is the mean velocity, T is the temperature,^ P is the 
pressure or stress tensor, and q is the heat flux. The trace of the pressure 
tensor defines the hydrostatic pressure p = jTr P, which in a dilute gas is 
related to the density and the temperature through the ideal gas equation 
of state p — fiT. In (2), > is the (local) "coohng" rate associated with 
the coUisional inelasticity and 7 > is the (local) "heating" rate associated 
with the external driving. If the collisions are elastic, the cooling rate vanishes 
(C = 0) and no energy needs to be externally injected into the system (7 = 0). 
On the other hand, in the case of inelastic collisions one has C > ^6 
will assume that the external source of energy compensates locally for the 
coUisional coohng, i.e., 7 = C- 

Now we particiilarizc to the gravity-driven Poiseuille flow. In the planar 
geometry, the gas is enclosed between two infinite parallel plates normal to 
the y-axis and the acceleration of gravity g = —gx is applied along a direction 
z parallel to the plates. A longitudinal section of the system is sketched in 
Fig. 1, where the variable s refers to the coordinate y. In the stationary laminar 
state the macroscopic variables depend on the variable y only, so the balance 
equations (1) and (2) yield 



dy ~ ' dy ~ '^^ ' dy ^ dy ~ ' ^ ' 



where p = m,n is the mass density and use has been made of the assumptions 
C = 7 = (elastic collisions) or C = 7 > (inelastic collisions). Equations 
(3) imply that the normal stress Pyy is uniform across the system. Note that 

* We have set the Boltzmann constant equal to 1, so that T (which is usually 
termed "granular temperature" in the case of inelastic collisions) has dimensions 
of an energy. 
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neither the other two normal stresses P^x and P^z nor the longitudinal com- 
ponent Qz of the heat flux appear explicitly in the balance equations (3) . 

The Poiseuille flow can also be defined in a cylindrical geometry. In that 
case, the gas is inside an infinitely long straight tube of circular cross section, 
parallel to the z-axis. Figure 1 also represents a longitudinal section of the 
sy stem, ex cept that now the variable s refers to the radial coordinate r = 
\/x^ + y^. Assuming again a stationary laminar flow, the only relevant spatial 
variable is r. Expressing the pressure tensor and the heat flux in cylindrical 
coordinates, the balance equations (1) and (2) reduce to [9, 11] 



— {rPrr)=P^^, -—{rPrz) = -pg, Pr.r—^ + — {rqr)=Q. (4) 



Our aim is to get the hydrodynamic profiles p(s) = n{s)T{s), Uz{s), and 
T{s) (where s = y and s = r for the planar and cylindrical geometries, respec- 
tively) in the hulk domain, i.e., around s = 0, by a perturbation expansion 
through second order in g. 

3 Navier— Stokes Theory 

Equations (3) (planar case) and (4) (cylindrical case) do not constitute a 
closed set of equations. In a purely hydrodynamic approach, they must be 
complemented by constitutive relations expressing the fluxes P and q in terms 
of the hydrodynamic fields and their gradients. 

If the gravity field were switched off (5 = 0), the physical solution to (3) 
and (4) would correspond to an isotropic and uniform state, i.e., = VT = 
Vu = 0, P = p I (where I is the unit tensor), and q = 0. This is not but 
the equilibrium state in the case of elastic collisions and a uniformly heated 
non-equilibrium steady state in the case of inelastic collisions. On the other 
hand, the existence of the constant field g, no matter how small it is, induces 
gradients along the directions normal to the boundaries. It seems then quite 
natural that a Navier-Stokes (NS) description might be appropriate, at least 
for asymptotically small values of g. As will be seen, this expectation turns 
out not to hold true. 

In a monodisperse gas, the NS equations follow from Newton's viscosity 
law and Fourier's heat conduction law. 



where the dagger denotes the transpose, 77 is the shear viscosity, k is the 
thermal conductivity, and /U is a coefficient that vanishes if the collisions are 
elastic but is in general non-zero in granular gases. The expressions of the 
transport coefficients rj and k in terms of the temperature and the mechanical 
properties of the particles are well known in the case of elastic systems [15]. 




(5) 
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Moreover, r], k, and jj. have been derived for inelastic hard spheres [16, 17] 
and inelastic Maxwell particles [18] in the free cooling case as well as under 
heating. 

Since in the Poiseuille flow the gradients are normal to the flow direction, 
it follows from (5) that the three normal stresses are identical. In Cartesian 
and cylindrical coordinates, 

Pxx ~ ^yy ~ Pzz ~ P ) -^<^<^ ~ ^rr ~ ^zz ~ P 5 (6) 

respectively. Equations (3) or (4) then imply that p = const, so that the 
(modified) Fourier law in (5) becomes 

q = -k'VT , (7) 

where we have called k' = k — n^/T. Moreover, there does not exist a longi- 
tudinal component of the heat flux. 

When the NS constitutive relations (5) are used in the balance equations 
(3) for the planar case, one gets the closed set of equations 

d du^ _ mpg d ,dT_ f duA^ 

Equations (8) give a parabolic-like velocity proflle, that is characteristic of the 
Poiseuille flow, and a quartic-like shape for the temperature proflle. Strictly 
speaking, these NS proflles are more complicated than just polynomials due 
to the temperature dependence of the transport coefficients. The knowledge 
of such a dependence is needed in order to get the NS hydrodynamic profiles 
from (8). However, if one is interested in those profiles to second order in g, 
the precise dependence of r] and k' on T is not needed and one can replace 
T] ^ rio and k' Kq.^ More explicitly, the solution to (8) through order g"^ is 

uz{y) = uo + ^y^ + 0{g'), p{y) = po , (9) 



2r?o" 



T{y) = To ( 1 - / ) + 0(/) . (10) 

The shear stress and the heat fiux are 

Pyz{y) = -Pogy U + j}^' y') + 0{g') , (ii) 

qy{y) = ^y' + 0{g'), qz{y) = 0, (12) 

* Henceforth the subscript attached to a quantity will denote the value of that 
quantity at the mid plane j/ = or at the symmetry axis r = 0. 
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respectively. Note that the knowledge of the temperature profile to 0{g'^) 
allows us to get Pj.y to 0{g^) through (3). Equation (10) shows that the NS 
solution predicts that the temperature has a maximum at y = 0. 
In the cylindrical case, insertion of (5) into (4) yields 

_i a / du,\ mpg _i d f , dT\ ( du,V 



dr \ dr J T ' dr \ dr J \ dr 

Again, these two equations can be easily solved to second order in g: 

Mr)=uo + ^r^ + 0{g^), p{r) = po , (14) 
r(r) = To f 1 - i- + 0(/) . (15) 



64 vok'oTo 

The corresponding expressions for the shear stress and the heat flux are 

^-(^)-fK^+ife;^^^)-^(^^)' (^«) 

Qrir) = ^r^ + 0{g'), g.(r) = , (17) 
16% 

respectively. Equation (15) implies that the temperature is maximum at r = 0. 



4 Kinetic Theory Description. Elastic Collisions 

In a kinetic theory approach the relevant information of the system is conveyed 

by the one-particle velocity distribution function /(r, v, t). Its temporal evolu- 
tion is governed by the non- linear Boltzmann equation [15], which in standard 
notation reads 

|/ + v.V/ + g-|^/ = J[/,/], (18) 

where J[f, f] is the collision operator for elastic collisions and we have as- 
sumed that no external forcing exists, apart from gravity. The influence of 
the interaction potential appears in J[f,f] through the dependence of the 
collision rate on the relative velocity and the scattering angle. The hydrody- 
namic variables and their fluxes are just velocity moments of the distribution 
function; 

{n, nu,p, P, q} = / dv {l, v, jV, mVV, f V^'v} / , (19) 

where the peculiar velocity V = v — u has been introduced as the velocity of 

a particle relative to the flow velocity. The exact balance equations (1) and 
(2), with (^ = 7 = 0, follow from (18) by taking velocity moments. 
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4.1 Planar Geometry 

In the case of the stationary Poiseuille flow in a channel, the Boltzmann 
equation (18) becomes 

Strictly speaking, (20) must be supplemented by the appropriate boundary 

conditions describing the interaction of the particles with the plates. However, 
we are interested in the bulk region of the system, sufficiently far from the 
boundaries. Therefore, we will look for the Hilbert-class or normal solution to 
(20), namely a solution where all the spatial dependence of the distribution 
function takes place through a functional dependence of / on the hydrody- 
namic fields n, u, and T. 

Since the hydrodynamic variables of the gas and the associated fluxes 
are the first few moments of the distribution function /, it is convenient to 
consider the hierarchy of moment equations stemming from the Boltzmann 
equation (20). A moment of an arbitrary degree = fci + A;2 + fcs is defined as 

MkMMa) = J dy^V,'^Vy'W,'^f{y,v;g) . (21) 

Because of the symmetry properties of the problem, Mki,k2-ka is a function of y 
of the same parity as ^2- Seen as a function of g, M^^^^zM has the same parity 
as fca. Finally, Mk^^k2,k3 = if fci = odd. Prom (20) one gets the hierarchy 

d ( duz \ 

^Affei,fe2+i,fe3 + ^3 ( -^Affci,fc2+i,fe3-i +5-^fei,fc2,fc3-i ) = JkiMM > (22) 

where JkiMM is a collisional moment. In general, it involves all the velocity 

moments of the distribution, including those of a higher degree, and its explicit 
expression in terms of those moments is unknown. An important exception is 
provided by the Maxwell interaction potential, in which case the collision rate 
is independent of the velocity and, as a consequence, JkiMM '^^^ ^le expressed 
as a bilinear combination of moments of the same or smaller degree. 

Even in the case of Maxwell molecules, the hierarchy (22) couples moments 
of a certain degree k to moments of degree fc + 1.^ However, the problem can be 
solved by a recursive scheme [6]. First, the moments are expanded in powers 
of g around a reference equilibrium state parameterized by uq, po, and Tq: 

oo oo 

w^(2/;5) = wo + X!"i^nyy > p{y;9)=Po + Y.P^^^'(y)9^ ^ (23) 

£=1 i=2 



' In what follows, wc will use the romari boldface k to denote the triad {/ji, /C2, kz} 
and the italic lightface k to denote the sum fci + fe2 + fcs. 
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Due to symmetry reasons, ui (y) = if ^ = even, p^^\y) = T^^\y) = if 
£ = odd, and M^^ (y) = if fes + f = odd. The second step consists of the 
ansatz that the coefficients have just a polynomial dependence on y, namely 



J=2 



(25) 



(26) 



(id) _ 







3=2 j=0 

where A^i = 1 and Ng = 2i for £ > 2. Symmetry implies that Mj 
if k2+j = odd. The numerical coefficients u^z'^^\ p(^.2j)^ T{e,2j)^ ^nd 
are determined recursively by inserting (23)-(26) into (22) and equating the 
coefficients of the same powers in g and y in both sides. This yields a hierarchy 
of linear equations for the unknowns. This rather cumbersome scheme has 
boon solved through order g^ in [6] . The results for the hydrodynamic proffles 
and the fluxes are 



u.{y)=uo + ^y' + 0{g') 



p{y) = Po 



T{y) = To 

Pzz{y) =Po 



1- 



12r]oKoTo 



y^ + Cr 



mg 



1 + 



yy 



Po 



Pyz{y) = -pogy 



9 9 



mg 

-c' 



y 



y^ + c. 



pmW 



pI 



Po 



0{g') 



+o(/) 

0{g'), 

^o(/), 



Ct ( TTT-g 



■0{g') 



1 1 

Po9 „,3 



<ly{y) = + 0{g^) , q,{y) = C.mgKo + 0{g') . 



(27) 

(28) 
(29) 
(30) 
(31) 
(32) 



The numerical values of the coefficients Cp, Ct, C^, C'^, and Cq are displayed 
in Tabic 1. Comparison of (27) (32) with the NS predictions (6) and (9) 
(12) shows that the NS approach fails to account for the non-zero values of 
those five coefficients. At a qualitative level, the main contrast appears for the 
temperature profile. While, according to the NS equations, the temperature 
has a maximum at the mid plane y = 0, the Boltzmann equation shows that. 
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because of the extra quadratic term headed by Ct, the temperature actually 
presents a local minimum Tq at j/ = surrounded by two symmetric maxima 
T^max at 2/ = ±2/max, where y^ax = \/^Crrm^/Po- The relative height of 
the maxima is (T^ax - Tb)/To = 2CT(52/max/wo)^, where vq = sfxT^jm is 
the thermal velocity at y = 0. If one defines an effective mean free path as 
A = 16v'2??Kr/157r/5p,'^ one can write y^ax = (15^/5^FCV/16)Ao ~ 3.74Ao, 
(T^a. - To)/To = (11257rC2/I28)(gAo/i;g)2 ~ 2%.m{g\^lvlf . 

Table 1. Numerical values of the coefficients Cp, Ct, C^, C^, and C,, both in the 
planar and the cylindrical geometries, according to the Navier-Stokes description 
(NS), the Boltzmann equation for Maxwell molecules (B-M), and the Bhatnagar- 
Gross-Krook kinetic model for ain- iutcractiou (BOK) 



planar cylindrical 



coefRcient NS 


B-M BGK 


NS 


B-M 


BGK 







1.2 1.2 





0.3 


0.3 


Ct 





1.0153 0.76 





0.19429 


0.14 







6.4777 13.12 





3.4776 


7.36 







6.2602 12.24 





1.7388 


3.68 


c. 





0.4 0.4 





0.4 


0.4 



Equations (27)-(32) are exact in the context of the Boltzmann equation 

for Maxwell molecules. For more realistic potentials (e.g., hard spheres) the 
hierarchy (22) cannot be solved recursively, so that the problem must be ad- 
dressed by means of approximations. One theoretical possibility is to consider 
the Chapman-Enskog expansion and retain as many terms as needed to get 
the profiles through order g"^, but this is not very practical. The coefficients 
Cp and Cq are already captured by the Burnett description [6, 7], and so 
the Maxwell values Cp = ^ = 1.2 and Cq = | = 0.4 can be expected to be 
good approximations for other potentials. On the other hand, the determi- 
nation of Ct, Cro, and C'^ requires the consideration of super-Burnett and 
super-super-Burnctt contributions [6]. A second approach is Grad's moment 
method [19]. The thirteen-moment approximation [3] yields (27)-(32) with 
Cp = |, Cq = |, Ct = 55 = 0.56, and C^ = C'^ = 0, irrespective of the 
interaction potential. On the other hand, the nineteen-moment approxima- 
tion [8] gives Cp = f29p ~ 1.214 and Ct = fffH 0.992 for hard spheres, 
while it predicts Cp = | and Ct = if = 1-04 for Maxwell molecules. Compar- 
ison with the exact results Cp = | and Ct = 1.0153 for Maxwell molecules 
suggests that the nineteen-moment approximation is rather accurate, although 
it overestimates Ct by a few percent. 

An alternative route to get reasonable estimates with relatively much less 
effort than in the case of the Boltzmann equation consists of using a ki- 

® This definition guarantees that A = (V27rncr^)^^ for hard spheres of diameter a. 
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netic model equation. The prototype kinetic model is the one proposed by 

Bhatnagar, Gross, and Krook (BGK) and, independently, by Welandcr and 
Kogan [20]. In this model, the Boltzmann collision operator is replaced by a 
relaxation-time term towards the local equilibrium state: 



The BGK model has proven to be quite reliable in non-Newtonian shear flow 

problems [21]. The solution to the BGK equation for the plane Poiscuillc 
flow has been explicitly obtained through order [5]. The results strongly 
suggest that the series expansion is only asymptotic, so that from a practical 
point of view one can focus on the first few terms. The results agree with the 
profiles (27)-(32), except that the numerical values of the coefficients Ct, C^, 
and C4j differ from those derived from the Boltzmann equation for Maxwell 
molecules, as shown in Table 1. The coefficient Ct in the BGK model is about 
25% smaller than in the Boltzmann equation, while and C'^ are about 
twice larger in the former than in the latter. Interestingly enough, the BGK 
value Ct = H = O-^S agrees quite well with Monte Carlo simulations of the 
Boltzmann equation for hard spheres [2], even though the nineteen- moment 
method predicts Ct < 0.99. 

4.2 Cylindrical Geometry 

In the cylindrical Poiseuille flow, the Boltzmann equation (18) becomes 



The algebra is now more involved than in the planar case, but still the hier- 
archy of moment equations can be recursively solved for Maxwell molecules 
through an expansion in powers oi g. To second order, the solution is [11] 




(33) 




(34) 





(36) 




(37) 




(38) 
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1.08 
1.04 
1.00 
^ 0.96 
0.92 
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■ ■ ■ NS ' ■ • . 
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-■- BGK 


A. 




B-M 


\ 








Fig. 2. Plot of (a) T{s)/To and (b) [P^^(s) ~ Pssis)]/po for gXo/vf, = 0.05 in the 
planar (s = y, thick lines) and cylindrical (s = r, thin lines) geometries, according 
to the NS equations {■■■), the BGK model (-•-•-), and the Boltzmann equation 
for Maxwell molecules ( — ) 



= 



1 + 



16r?o 



192r]oKoTo 



4 1 

r H — — 



mg 
To 



+ 0(9"), (39) 



qz{r) = CgmgKo + 0{g^) 



(40) 



The numerical values of the coefficients Cp, Ct, C^, C'^, and Cq for this 
case are again displayed in Table 1. As happened in the planar case, the 
temperature exhibits a non-monotonic profile: it has a local minimum Tq 



at r 



and its maximum value T, 



is reached at r = rmax: where 
(5V157rCT/4)Ao ~ 3.78Ao. The relative value of 



?*max = V32Ct?7oKoIo/po 

the maximum is (Tmax-To)/To = 2CT(5r-max/«o)^ = iS75TrC^/8)igXo/v^ 
5.56(.gAo/uo)^. It is noteworthy that, expressed in the same units, r 



,2\2 



is 



practically equal to j/max, while Tir 



To is about five times smaller in the 



cylindrical geometry than in the planar geometry. 

The BGK model has been solved for the cylindrical Poiseuille flow through 
order [9], suggesting again an asymptotic character of the expansion. The 
resulting profiles agree with (35)-(40), but with different values of Ct, C^, 
and C4j (see Table 1). It is interesting to note that the ratios between the 
BGK and the B-M values for Ct, Ct^, and C'^ are almost the same in the 
planar and the cylindrical geometries. 

Figure 2 shows the profiles for the temperature and for the normal stress 
difference Pzz — Pss (where Pss = Pyy and Pss = Prr in the planar and 
cylindrical geometries, respectively) for gXo/vQ = 0.05. 



5 Kinetic Theory Description. Inelastic Collisions 

Now we consider the plane Poiseuille flow and assume that the particles are 

smooth inelastic hard spheres of diameter a and coefficient of restitution 
a < 1. In the dilute regime, the cooling rate is approximately given by [16, 22] 
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C(«) = v^{l-a^), -^na^ {TrT/mf^ . (41) 

It is obvious that all the mathematical difficulties embodied in the Boltzmann 
equation for (elastic) hard spheres are further increased in the inelastic case. 
Therefore, in order to get explicit results with a reasonable amount of effort, 
it seems convenient to extend the BGK model (33) to the realm of inelastic 
collisions. The simplest possibility is perhaps [22, 23] 

J[fJ] - -PiaMf - /,) + . (V/) , (42) 

where l3{a) is a dimensionless function of the coefficient of restitution that can 
be freely chosen to optimize the agreement with the Boltzmann description. 
A simple choice is P{a) = ^(1 + a) [23]. Here, however, we will take (3{a) = 
|(1 + a) (2 + a), which makes the BGK shear viscosity agree well with the 
Boltzmann one [22, 23]. Good quantitative agreement between the kinetic 
model (42) and the Boltzmann equation has been found for the simple shear 
and Couette flows [24]. 

As said in Sect. 2, an external energy input is needed to compensate for the 
coUisional dissipation in order to reach a nonequilibrium steady state. In real 
experiments this is usually achieved by means of boundary vibrations of small 
amplitude and high frequency. However, this type of realistic heating through 
the boundaries is difficult to deal with at a theoretical level due to unavoidable 
boundary effects. These difficulties are overcome by assuming a bulk heating 
mechanism acting on all the particles simultaneously. The most commonly 
used type of bulk driving for inelastic particles consists of a stochastic force in 
the form of Gaussian white noise, which appears in the Boltzmann equation 
under the form of a diffusion term in velocity space [25]. In summary, our 
kinetic equation to describe the stationary plane Poiseuille flow reads [12] 

where the heating rate 7 is taken equal to C, as discussed in Sect. 2. 

The perturbative solution of (43) is seen to agree with the structure of 
(27)-(32), except that now the NS transport coefficients ry and k' = K — n^/T, 
as well as Ct, C-^, and C^, are functions of the coefficient of restitution. 
As said above, the choice /? ~ i(l + a)(2 + a) guarantees that the BGK 
shear viscosity 77 agrees with the one derived from the Boltzmann equation 
for inelastic hard spheres in the simplest Sonine approximation. However, as 
already happens in the elastic case, the (modified) thermal conductivity k' 
differs in both kinetic equations. In order to circumvent this problem, once 
the BGK solution is expressed in terms of rj and k' , we will take for those 
coefficients the expressions obtained from the Boltzmann equation. In the case 
of a granular gas heated by the white noise forcing, the transport coefficients 
are approximately given by [17] 
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a a 

Fig. 3. Plot of (a) (Tmax — To) /To for gXo/vo = 0.05 and (b) 2/max/Ao as functions 
of the coefficient of restitution a 



, = P 1 «'=i^ ^1 (44) 

While rj monotonically increases with inelasticity, k' starts decreasing with 
increasing inelasticity, reaches a minimum value around a ~ 0.4, and then 
slightly increases for a < 0.4. The expressions for Ct{o!), C^{a), and C!^{a) 
are 

^ 1 38 + 43r + 17r^ ^ M^-a') , . 

"^^^"^ 25 (1 + C*)(2 + C*) ' ^ ~ /3(a) + ^(l-a2) ' ^ ' 

16 82 + 67C+8C! 12 102 + 87C + 13C" 

"^^"^ 25(1 + C)(2 + C*)^ ' 25 (1 + C*)(2 + C)' ■ 

Of course, in the clastic limit {a 1) the above coefficients reduce to the BGK 
values listed in Table 1. The coefficient Ct(q;) monotonically decreases with 
increasing inelasticity. However, the non-monotonic behavior of K'{a) makes 
the relative location of the maximum j/max/Ao, where Aq = (\/27moO"^)~^, as 
well as the relative height (Tmax — Tq)/Tq for a given value of gXo/vQ, exhibit 
a non-monotonic behavior. This is illustrated in Fig. 3. 



6 Conclusions 



In this paper wc have reviewed the main results derived within a kinetic the- 
ory framework for the steady Poiseuille flow driven by a uniform longitudinal 
force. The results are restricted to the bulk domain and are obtained by a per- 
turbation expansion through second order. We have considered both elastic 
and inelastic particles. In the first case, the planar and the cylindrical geome- 
tries have been considered and the results are presented from the Boltzmann 
equation for Maxwell molecules as well as from the BGK kinetic model for 
any interaction. In the case of granular gases (inelastic collisions), however, 
only the planar geometry and a BGK-like kinetic model have been addressed. 
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In all the cases, the profiles of the hydrodynamic fields and their fluxes 
exhibit important deviations from the NS predictions: the hydrostatic pressure 
is not uniform, normal stress differences are present, a component of the heat 
flux normal to the thermal gradient exists, and the temperature profile has a 
non-monotonic shape. In general, those deviations arc more important in the 
planar geometry than in the cylindrical one. For instance, the relative height 
of the maximum temperature, (Tmax — 7o)/7o, is about five times larger in the 
former case than in the latter, even though the location of that maximum is 
similar in both cases, i.e., ymax — ?'max- It is worthwhile noting that the BGK 
model succeeds in capturing the functional form of the profiles derived from 
a more fundamental Boltzmann description, albeit with some changes in the 
numerical factors. Thus, the BGK model underestimates the maximum value 
T'max, as well as its location ?/max or rmax- Concerning the case of inelastic 
collisions, the results show a weak influence of the coefficient of restitution a. 
In any case, for small and moderate inelasticities (say a > 0.5) there is a slight 
decrease in the quantitative deviations from the NS profiles as inelasticity 
grows, while the opposite behavior takes place for high inelasticity (a < 0.5). 
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